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Abstract 

The equations describing a two-component cosmological fluid with linearized density perturba- 
tions are investigated in the small wavelength or large k limit. The equations are formulated to 
include a baryonic component, as well as either a hot dark matter (HDM) or cold dark matter 
(CDM) component. Previous work done on such a system in static spacetime is extended to reveal 
■ some interesting physical properties, such as the Jeans wavenumber of the mixture, and resonant 

mode amplitudes. A WKB technique is then developed to study the expanding universe equations 

> 

in detail, and to see whether such physical properties are also of relevance in this more realistic 

scenario. The Jeans wavenumber of the mixture is re-interpreted for the case of an expanding 

background spacetime. The various modes are obtained to leading order, and the amplitudes of 

the modes are examined in detail to compare to the resonances observed in the static spacetime 

Q_i' results. It is found that some conclusions made in the literature about static spacetime results 

2 ■ cannot be carried over to an expanding cosmology. 
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I. INTRODUCTION 



The analysis of cosmological perturbations in the Newtonian limit is a well studied prob- 
lem in theories of structure formation, and it may be supposed that there is little left to 
learn from this theory. Most of the effort has gone into the study of the one-component 
cosmological fluid equations, and the results have been well expounded in many standard 
texts . There is, however, still a wealth of problems remaining in the detailed 

analysis of two-component cosmological fluids and their linearized gravitational perturbation 
modes. In particular, if pressure effects are included so that the Jeans instability becomes an 
issue, the equations present a considerable analytic challenge, and a range of new physical 
effects become apparent. Some of these effects have been studied in the contrived case of 
a static spacetime background ((J In this scenario there is no expansion, so that the 
mathematics is considerably simplified, and solutions can easily be found. This is useful to 
gain some qualitative idea about physical phenomena observable, but to gain a true picture 
in a cosmological context, the expanding background spacetime given by the Friedmann- 
Robert son- Walker cosmologies is required. 

There have been a variety of studies of the mult i- component cosmological fluid equations, 



a, 



ranging from some relatively specific applications under certain cosmological scenarios 
to a broad mathematical study and classification [lOj. A discussion of the application and 
validity of some of the equations mentioned in these previous studies, together with the so- 
lution of an unsolved set of two-component post-recombination equations, has recently been 
undertaken by the authors jllj . The system of equations described the interaction between 
a dark matter and baryonic component in the Newtonian regime (density fluctuations on 
scales well within the Hubble radius). A series expansion of the solutions for small wavenum- 
ber k (large scales) was presented. This allowed comparison with some of the previous work, 
in particular with the Meijer G-function classifications given by jioj . This region of fc-space 
is also interesting because it is the region in which the Jeans instability is known to occur. 

In this paper we wish to complete this study by examining the large k asymptotic region of 
the solutions. Such a study is worthwhile, in order to make contact with the static spacetime 
results of Although not realistic as cosmological solutions, these results displayed a 
number of little known physical phenomena associated with the linearized modes, which we 
wish to expand on here. The techniques required to analyze the expanding universe solutions 
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are also of interest in their own right mathematically, where a generalized WKB method 
will be expounded. It is possible to make a comparison with the work done in cosmological 



plasma physics in an Einstein-deSitter background [13j,|lJ|. This is interesting because of 
the mathematically very similar form of fluid equations for both type of systems, which is 
due to the similarity of the electromagnetic and gravitational forces. Thus mathematical 
techniques employed in the analysis of plasma equations will be useful in this paper, and 
give clues as to how to proceed with some challenging mathematical analysis of gravitational 
density perturbation modes. 

The paper is to be organized as follows. The relevant equations will be introduced in 
Section 2. The discussion will then be focused in Section 3, by reconsidering the two- 
component modes in a static spacetime. This investigation is by necessity of a qualitative 
nature, but gives a useful introduction to the concepts and interesting physical effects not 
found in the standard one-component analysis. The work of [3] will also be extended. The 
expanding universe baryonic and dark matter equations will then be considered in Section 4. 
The short wavelength (WKB) approximation will be utilized to complete the study of these 
equations initiated in The relevance of the previous work on static spacetime systems 
is revealed through this analysis. This will allow meaningful conclusions to be drawn about 
this whole area of study, and point to where promising future work may lie. These aspects 
are discussed in Section 5. 



II. THE GOVERNING EQUATIONS 

A broad survey of the Newtonian cosmological perturbation equations under various 
cosmological scenarios was given in [llj]. This paper also gave a detailed derivation of the 
equations of interest for present analysis. We will not repeat such a detailed discussion here, 
but directly introduce the relevant equations. 

The starting point is the equations for an n-component system of nonrelativistic species, 
as derived in all the standard texts. Given a density perturbation 5, of the i-th component 
of the mass density pf. 

= (2-1) 
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it may be decomposed into its Fourier plane wave modes with wave vector k 

Si(r,t) = -ip J 5 kl (t)ex P (-tk-r)d 3 r. (2.2) 

Here r is the physical spatial coordinate, and t is cosmic time. Using the Eulerian equations 
of motion describing a perfect fluid, a set of coupled second order equations for the Fourier 
modes 5i(t) (where we now drop the subscript k) are achieved: 

^ + 2 ^ + ^ i = 4?rG ^ M ' l = l ^--- n - ( 2 - 3 ) 

i=l 

Overdots will denote derivatives with respect to t. 

The above equations contain the universe expansion factor a and sound velocities Wj. 
We use the expression "sound velocity" fairly loosely. The parameters Vi could also denote 
a general velocity dispersion for a collisionless fluid. To be able to solve the equations, 
the time dependence of the physical variables needs to be made explicit. We will adopt 
the convention that barred variables will denote comoving quantities, independent of time. 
Thus the definition of comoving wave number k = ak arises naturally. An assumption must 
also be made about the scaling of the sound velocities i>j. In the post-recombination era, 
the adiabatic speed of sound follows the behavior v oc a -1 , so we will introduce the time 
independent quantities = avi. The total background energy density p Q can also be made 
independent of time by the definition p = a 3 po- This enables us to introduce the useful 
parameter = Pi/po, the fraction of mass density contributed by species i. 

We will not go into detail on how the equations are transformed into their simplest 
form here (instead see [11]), but briefly describe the important points. The equations are 
first transformed so that a is the only explicit temporal variable. This allows parameters 
specifying particular large scale cosmological dynamics to enter the equations, namely the 
cosmological constant A and spatial curvature k c . In this paper, a background Einstein 
deSitter cosmology will be employed (k c = 0, A = 0). Although this model has been ruled 
out with high confidence by current observations, it is sufficient for our purposes. We wish 
to study some important physical process without additional complications. 

It is found that the following important parameters arise: 

- 2 _ 4nGp - 2 _ AnGpo 

B - —o ' K D - ~2 • 

U B U D 

They resemble the comoving Jeans wavenumbers for each component taken separately. The 
strict Jeans one-component wavenumbers are given by replacing p with p^ in (12. 4|) (see 



next section). The wavelength parameters defined above indicate whether gravity ki > k 
or pressure support k{ < k dominate the dynamics, and thus whether the region of fc-space 
under consideration is Jeans unstable. It may also be noted that the relation e B + = 1 
holds. 

It is useful to define the parameters 

K B = A K D = A, (2.5) 
k B k D 

for a clear dimensionless partitioning of parameter space. < 1 corresponds to the Jeans 
unstable region in the single component analog of the equations, and Ki > 1 to the acoustic 
region. The cosmological fluid equations are finally written in terms of the variable \ = 
a -1 / 2 , to give the canonical form of the system of differential equations to be studied in this 
paper: 

5' B + 6(k b - 5§) 5 B - 6 -^6 D = 0, (2.6) 
S' D + 6 (k% - 8 D - 6 -^5 B = 0. (2.7) 

The prime denotes differentiation with respect to \- 

These equations bear a strong resemblance to the equations of an electron-proton cos- 
mological plasma studied in Q| [equations (4.8) and (4.9) of that paper]. Considering the 
mathematical similarity between the electromagnetic and gravitational forces, this was to 
be expected. A manifestation of this fact is the close resemblance between the dispersion 
relation for the sim ple one-component Jeans instability and Langmuir modes. The tech- 
niques employed in will be adopted and developed further to the current problem. In 
particular, some general WKB techniques will be extended. This will also indicate further 
results obtainable in cosmological plasma physics. 

We now digress to an analysis of the static spacetime perturbation equations to introduce 
some new physical phenomena, which are to be scrutinized for their applicability in an 
expanding universe. 
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III. THE STATIC TWO-COMPONENT PROBLEM 



A. Eigenvalues and Eigenvectors 

The static spacetime results for a two-component fluid are well understood, though receive 
little attention in standard linearized structure formation theory, which aims to produce the 
power spectrum of density perturbations. We will extend the current results to facilitate 
understanding the general expanding universe scenario later. This section aims to develop 
some concepts in a relatively simple setting. Previous work on the static problem has been 
done in (Jill. We will in particular rely quite heavily on the notation and results of []| 
in this section. Although the solutions are unrealistic as an application to cosmology, they 
display some similar qualitative features, and allow an exposition of the basic physical ideas 
without the complication of spacetime expansion being introduced. The static nature of 
the spacetime simplifies the mathematics greatly, and is thus useful in understanding the 
general problem. 

In this section, there is no need to refer to barred (comoving) physical quantities, and 
all physical variables may be assumed to be constant in time, unless otherwise specified. 
With the expansion parameter a set equal to unity, the general fluid equations ()2.3|) may be 
written as 

5d + iy 2 D k 2 - W D )5 D - W B 5 B = 0, (3.1) 
6 B + (v 2 B k 2 - W B )5 B - W D 5 D = 0, (3.2) 

with an overdot denoting differentiation with respect to t, and Wi = AnGpi. A study of 
the behavior of the solutions to these equations is most readily undertaken by reducing the 
system to a first order autonomous dynamical system, undertaken in Q. To analyse the 
dynamical system, a solution needs to be found for the state vector 

x = (xx, x 2 , x 3 , x 4 ) T = (5 D , 5 D , S B , S B ) T - (3.3) 

We just state the results here. 

The most important feature discovered in [?| was the existence of a parameter dependent 
critical point of the dynamical system given by 

k = k M = k B + k D = g I ~ ' (3-4) 

V B V D 



where k B and kn have been defined in terms of the density and velocity parameters of 
each matter component, and are slightly different from k B and hp defined in (12. 4 j) . The 
special value of the wavenumber ku-, may be thought of as the Jeans wavenumber of a two- 
component fluid (the mixture wavenumber). It comprises the Jeans wavenumbers of each 
fluid taken separately, but it is possible to show that k = ku is the only physical quantity 
which indicates an instability — both k = k B and k = k B have no such interpretation for the 
coupled two-component case. 
With solutions of the form 

4 

x(i) = y]a i exp(Ait)£ i , (3.5) 

i=i 

where the aij are amplitude functions dependent on k and determined by initial conditions, 
the solutions for the eigenvalues Aj and eigenvectors ^ of the dynamical system are respec- 
tively 

(\l = -a 2 = ^f+^/pTfg 
\ A 3 = -A 4 = j=Jf- \[f r +fg 

with 

/(A;) = w B + W D -k 2 (v 2 B + v 2 D ), 
g{k) = k 2 (W B v 2 D + W D v 2 B )-k 2 v 2 B v 2 D , 

and 

& = (AA i ,A,A i ,l) r , 2 = 1,2,3,4, 

with 

U = ft = sfe (HTFTlp^) 

and 

fc(fc) = W D -^ B + A; 2 (t;|-4). 
For calculational purposes, we note that 

/i 2 + 4W B W D = / 2 + 4c/. 

An examination of the real and imaginary parts of the Aj will show that the Ai and A2 modes 
represent acoustic oscillations for k > ku and growing and decaying modes for k < ku- The 
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(3.6) 

(3.7) 
(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 



A3 and A4 modes however always represent acoustic oscillations. The exponential nature of 
the solutions indicate that the growing and decaying modes do not have the typical power law 
behavior exhibited by expanding universe solutions, however the solutions exhibit the correct 
qualitative behavior in the regions below and above the critical point given by k = ku- 

To gain a feel for the properties of the above eigenvalues and eigenvectors, which allows 
us to make direct contact with the physics of the solutions, we study the quantities in 
various asymptotic regimes, and examine some plots. This will be effectively facilitated if 
the quantities are reparameterized in terms of some dimensionless variables. We need only 
consider the eigenvalues Ai and A3. To indicate the nature of the dark matter, the sound 
velocities may be coalesced into the single variable 



2 



V 2 = -f. (3.13) 



v 



Then V <C 1 corresponds to HDM while V 3> 1 corresponds to CDM. We also introduce the 
quantities en and e B as used elsewhere in the paper. In this context, they may be defined 

as 

W D W B 

€D = w B + w D > €B = w B + w D (3 - 14) 

We also parameterize the wavenumber dependence in units of the mixed Jeans wavenumber; 
thus we define 

K M = ^- (3.15) 

% 

It then follows that the eigenvalues may be written as 

Ai,3 = j=(W B + W D )W [l-Kl- (e D V 2 + K 



V2 
± 

+ 



2 

a; 



4(y + eDV) 2 Kh(l-Kh) 



1/2 



(3.16) 
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These expressions may be expanded for small and large Km- The results are: 
Ai ~ (W B + W D ) 1/2 



1 + ^(-£b + 4-^ 



1 

2e B e D V 2 - e D V* + e 2 D V*)K 2 M + ■■■], K M < 1 (3.17) 



~ %K M ^={W B + Wd) 1 ' 2 {1 + e B /V 2 + e D V 2 



- [(1 + e B /V 2 + e D V 2 ) 2 - A(e B /V + e D V) 2 } , K M » 1, (3.18) 

A 3 ~ iK M (W B + W D ) 1/2 (e B /V 2 + e D V 2 ), K M < 1 (3.19) 

~ iK M -^(W B + W D ) 1/2 {1 + e B /V 2 + e D V 2 
v2 

+ [(l + e B /V 2 + e D V 2 ) 2 -4(e B /V + e D V) 2 ] 1/2 y /2 , K M > 1. (3.20) 

These expansions confirm the earlier statement, whereby the A3 (and equivalently A4) modes 
display acoustic oscillations at all wavelengths, whereas the Ai (and equivalently A2) modes 
undergo a Jeans instability to growing (decaying) modes for K M < 1. It is also evident that 
at very large wavenumbers (small scales) the acoustic oscillations have a very large frequency, 
growing in proportion to the wavenumber, whereas for very low wavenumbers the A3 modes 
behave in a very slowly varying oscillatory manner, the frequency again being proportional 
to the wavenumber. In this regime the Ai and A2 modes comprise exponentially growing or 
decaying perturbations over an almost wavenumber independent timescale, approximately 
equal to (W B + Wd) 1 ^ 2 - These properties are illustrated in FiglH where the absolute values 
of the eigenvalues are plotted as a function of Km for a variety of H/CDM scenarios. The 
values e B = 0.1 and sd — 0.9 have been used in the plots, which is a fairly typical proportion 
of baryonic and dark matter mass density expected in the Universe. 



B. Interesting Scales 

It is evident that at certain scales the eigenvalues undergo some qualitatively interesting 
changes, which have been marked on the plots by some arrows. For wavenumbers around 
K M = 1, the Ai eigenvalue drops very quickly to zero, indicating the Jeans instability, but 
the A 3 eigenvalue displays uniform behavior in this region. There is another interesting scale 
in the Km < 1 region for small V. The physical motivation for this scale was discussed in 
0. It corresponds to a critical wavenumber kc, defined to be when the frequencies of each 
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component taken separately coincide, i.e. when 



v B k l -W B = v l D k l - W D . 



(3.21) 



The wavenumber k = kc is consequently given by 



k 



c 



W D -W B 

V 2 D -V 2 B 



1/2 



(3.22) 



It is interesting that the importance of this scale is only apparent for small V, where a 
sudden increase in the magnitude of Ai and decrease in the magnitude of A3 is apparent. 
For all V > 1, the plots would be almost identical to the displayed plots of V = 1000 in 
Figffl 

To gain a better understanding of this behavior, it is useful to convert kc into units of 
fcjw, which are the plotting units of all the figures. Thus 



Kmc = 



c 



1/2 



(3.23) 



k M V e D ~ £b ~ £dV 2 + e B /V 2 
It is interesting to compare this quantity to the individual Jeans instability scales for each 
fluid taken separately: 

k r < - ^ 1/2 



K M d = 
Kmb = 



ku 
k D 
ku 



e B /V 2 + e D 



1/2 



(3.24) 
(3.25) 



^B + e D V\ 

A better qualitative feel for these scales is facilitated by considering their expansions in the 
HDM and CDM regimes. For HDM, with V < 1 we find 



K 2 



K 2 



K 2 



V ' z [ — -1 



_ 1 ) v 2 + 0(V 4 

^B 



v 2 ^ 

t-B 



(-D- 



l-—V 2 + 0{V A ) 



(3.26) 

(3.27) 
(3.28) 



From this it may be concluded that if 
1. e D > e B , then K MC ~ K MD , 



2- e D > e B , then K MC < 1, 



3. e B > e B , then K MC is imaginary (no physical significance). 
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This last point is also borne out by the original definition (|3.22|) . where it is seen that for 
kc to be real, the dominant component must also be the hotter component. For CDM, with 
V ^> 1 the corresponding relations are given by 



K 2 - 


v~ 


2 ( £b 






K 2 - 


1 - 








CD 


K 2 - 


v~ 


2 e B 






CD 



CD 



1 v- 2 + o(v~ 



r-4\ 



(3.29) 
(3.30) 
(3.31) 



i-^y- 2 + o(K- 4 ) 

This shows that if 

1. e B > e D , then K MC ~ K MB , 

2. e B > e D , then K MC < 1, 

3. 6d > e_B, then Kmc is imaginary (no physical significance). 

The position of the arrows in Fig^bear out the above relations, as do the arrows in Fig|21 
and FigOlto be discussed more below. 

In conclusion, the fact that all eigenvalues have been plotted for the values e B = 0.1 and 
tr> = 0.9 means that the scale kc is only physically relevant for HDM. This is why all plots 
for V > 1 (CDM) are so similar. Given that the real Universe is now considered almost 
certainly CDM dominated, it is doubtful as to whether this potentially interesting physical 
effect given by the equations has any discernible effect on structure formation scenarios in 
pure CDM models. A H+CDM model may give similar interesting results as discussed here. 
This however is a three-component problem beyond the scope of this paper, but may be 
considered a worthwhile topic of research for future work in this area. 



C. Other Qualitative Behavior 

We briefly consider some other features of the solutions, to help better understand the 
mathematical properties. The qualitative features we wish to explore are the same for both 
the Ai and A3 modes, so we will concentrate only on the Ai modes here. The eigenvalue Ai 
is plotted for different values of e B and in Fig |2] and Fig|21 An interesting feature of these 
figures is that there is a one-to-one correspondence between each of the four plots in one 
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figure to a particular plot in the other figure, yet each plot corresponds to different physical 
parameters in each figure. This property highlights the symmetry of the eigenvalues. If the 
original analytic expression (J3.16|) for Ai is examined, it is clear that the expression retains 
an identical form if and e# are interchanged together with V and 1/V. Real values 
of Kmc have also been marked in. They indicate when the critical scale kc is physically 
relevant. Related to this property is the fact that in FigOall plots for €d < 0.1 and e# > 0.9 
are almost identical to the values of €d = 0.1 and e# = 0.9. Analogously, in FigJ3]the same 
may be said for all plots > 0.9 and e# < 0.1. 

We now turn to study the behavior of (3\ and fa, which give an indication of the relative 
proportion of baryonic and dark matter in each of the modes [see for example the X2 and X4 
components of the eigenvectors in Eq. (J3.9j) ]. In dimensionless variables, f3\ and (3% may be 
written as 



It is obvious from inspection that for all wavenumbers (3\ and (3^ are real valued and (3\ > 0, 
/?3 < 0. The almost symmetrical nature of the quantities are well illustrated in Fig0]for the 
opposing cases of H/CDM. Of interest here again is the scale Kmc, around which all the 
undergo an abrupt change. FigE] shows that either f3i or /5 3 will dominate very rapidly for 
increasing wavenumber, depending on the value of V. Stated more specifically, for the HDM 
scenario (V 1), baryons will dominate the Ai modes and dark matter will dominate the 
A3 modes, and vice versa for the CDM scenario (V ^> 1). We refrain from examining the 
asymptotics of (3\ and f3^ now, and leave that task to when we derive analogous expressions 
in the expanding universe scenario. The asymptotics will confirm the present qualitative 
discussion. 

D. Initial Conditions and Amplitudes 

Having completed a study of the behavior of the general solutions, we now turn to consider 
the effect of initial conditions. The study of the amplitudes of the various modes was studied 




(3.32) 
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with particular interest in []|, where the presence of a resonance was discovered at the scale 
kc- We extend that work here to consider the amplitude functions at a wider range of scales, 
and for CDM as well. This will be of relevance when the expanding Universe scenario is 
analyzed in the next section, where the Km > 1 range of scales needs to be considered. It is 
also of course of relevance from the fact that the Universe is believed to be CDM dominated. 

The amplitudes are fc-dependent functions, which also depend on the various constants 
in the problem. Consider the initial conditions at some time to given by the constants 
Xi{to) — x m- Here and henceforth, any variable subscripted with a (possibly together with 
other subscripts) denotes that quantity evaluated at t — to- As a reasonable simplifying 
assumption, the perturbations are assumed to start from rest, so that X\q = x^q = 0. The 
matter density perturbations may then be written in the form 

5d(t) = £ 2 o[Ci(e AlT + e- AlT ) + C2 cos(iA 3 r)], (3.33) 
8 B {r) = x 40 [C 3 (e AlT + e- AlT ) + U cos(iA 3 r)], (3.34) 

where the time origin has been shifted by the definition r = t — to- The amplitude functions 
Q are constructed out of the mode eigenvector functions fa, as well as the ratio of initial 
densities Q = x 40 /x 2 o- They are 

C3 " 2^^' C4 " 13^W (3 ' 36) 

Of particular interest is the fact that some of the Q display a resonance around the 
scale kc- This scale has of course previously shown its significance in the behavior of the 
eigenvalues and the fa. The fact that kc defines the scale at which the collapse times of 
the components taken separately coincide indicates that a resonance may well be expected 
to occur at this scale. The behavior of the various amplitudes over a wide range of scales, 
and in both the HDM and CDM scenarios are illustrated in FigslHHT^l The distinguishing 
feature of all the plots of HDM amplitudes is the rapid change of the functions around the 
scale kc- 

The a^ytic propel of Cl and to particular Cs are dl8 eu S8 ed at so,„e .engta in Q. 
Under some restrictive conditions (only HDM and certain initial values of Qo) it was shown 
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that (i would not obtain a resonance, whereas £3 would for 




) 



(3.37) 



In contrast no resonances are observed in the CDM scenario, but the amplitudes still undergo 



a rapid change around the scale K MB . Both K M = K MB and K M = K MD are always less 
than Km = 1, so that no significant behavior occurs in the long wavelength limit, a fact of 
some importance in later work. 

IV. THE SHORT WAVELENGTH APPROXIMATION IN THE EXPANDING 
UNIVERSE 

A. Matrix Formulation 

We are now prepared to tackle the most general equations formulated for the current 
problem, given by ()2.6|) and ()2.7|) . The solution to this system of equations cannot be 
classified by known analytic functions, so approximation schemes need to be implemented. 
This paper investigates a short wavelength approximation, which would be expected to probe 
the acoustic regime of the modes. A WKB type method may be employed for this. This 
is interesting, because through the derivation the explicit physical approximations required 
and type of solutions obtainable will naturally arise as a consequence of the method. A WKB 
method for coupled systems of equations in a cosmological plasma setting was expounded 
in We will further develop that method here for the current system, which is more 

complicated than anything considered previously. 

To begin with, the equations must be reduced to a first order system. Thus as in the 
static case ()3.3|) . we define 



x = (xi, x 2 , x 3 , x 4 ) T = (5' D , S D , 5 B , 8 B ) T . 



(4.1) 



The system may be written in matrix form 



x = 



Tx, 



(4.2) 
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with the definition 



-6Kl + ^ 

X 



6e 



B 





6e D 







X 




x l 





-6Kl 



6e B 



X 1 



(4.3) 



2 " " ^ n 

1 

The idea behind the method is to attempt to remove the coupling between equations as 
much as possible, hopefully relegating it to some lower order, which can then be dealt with 
by a suitable approximation. To this end we define the new matrices A and f such that 



x = Af. 



(4.4) 



A is chosen appropriately in order to diagonalise T. Then (14. 2 j) may formally be written as 



f = A TAf — A Af, detA^O. (4.5) 

To diagonalise T, we must first find its eigenvalues and eigenvectors. The structure of T 
is very similar to its static spacetime counterpart, so the four eigenvalues also have the form 
given by (|3.6|) . In the present case however, / and g are functions of x an d k, defined as 

f(x,k) = - 2 -Q(Kl + K 2 D ), (4.6) 
X 



3fi 

g( X , k) = -(Kle B + K%e D ) - 3QK 2 B K 2 D . (4.7) 

X ' 



Note that unlike the static spacetime results, where / contained an expression of the form 
Wb + Wd, no analogous expression exists here due to the physical constraint of the Einstein- 
deSitter universe, e# + eo = 1- The eigenvectors ^ corresponding the eigenvalues A« are also 
identical in structure to their static spacetime counterparts In this case the functions 

(5i are given by 



A = *=5i£( k -V^F^ (4 ' 9) 



with 

6 

~2 



h( X , k) = -{t D - e B ) + 6(K 2 B - K 2 D ). (4.10) 



X 
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It is worthwhile to point out here for the sake of calculations that 

,3Qe B e D 



S= a h 2 + A- 



X' 1 



(4.11) 



The eigenvectors may be used to form the diagonalising matrix 

A=(€x,* 2 ,*3>*4), (4-12) 
whose inverse exists. This enables the formal equation ()4.5|) to be written explicitly: 



(4.13) 



We have introduced eight new parameters here, all of which can be written in terms of Ai, 
A3, Pi and p 3 . In what follows, let us use the shorthand 

;2 




















' h 




ei 
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6e D ' 



■S 



+ 2x 2 (e D - e B ){Kl - iff,) + X *(K B - Klf. 



(4.14) 



The new parameters are defined by 

£3 
£5 



Pis 2A X ' 

a; 



2Ai 

2/3i3 V A! 

-Afi + ^ 

2/3i 3 V A 3 
A3 



2/3i 3 
/3i3 2A 3 : 

A 
2A3' 



(4.15) 
(4.16) 
(4.17) 
(4.18) 
(4.19) 
(4.20) 
(4.21) 
(4.22) 
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Let us consider the meaning of (|4.13|) more closely. Using the terminology of WKB theory, 
the first matrix on the right hand side will give us the leading order "control factor" — the 
fastest varying part of the solution, typically an exponential factor. This factor may indicate 
rapid oscillations for imaginary Aj, or rapid growth or decay for real Aj. The second matrix 
contains a collection of parameters, which determine further slowly varying behavior. For 
this to be true, the condition e. t <C Xj, Vz, j must hold. It is then possible to show that the 
four equations all decouple to leading order, and WKB solutions may be written down. The 
proof of this is quite involved, but it is worthwhile to pursue. A bonus of the proof is that 
through a careful consideration of the approximations required, some instructive physics is 
learnt along the way. 



B. The WKB Approximation Criteria 

In an attempt to decouple the equations, we need to consider more carefully the various 
criteria which constitute the condition ej <C Xj, Vi, j, which allow the WKB method to 
work. To begin with, we will assume the Aj's are of about the same order. Although their 
magnitude varies greatly over various scales, an examination of FigG] shows this assumption 
to hold fairly well in general. When we find, through the reasoning which follows, the precise 
scales of interest for a WKB approximation, we will see that this assumption is justified post 
facto. The essence of the WKB approximation is to assume 



< |Ai|, (4.23) 



that is, the eigenvalues vary slowly over the timescale which they define. In the current 
context, this corresponds physically to many oscillations in a universe expansion time for 
the acoustic region of /c-space, or a far shorter collapse time than a universe expansion time 
for perturbations in the region where modes are Jeans unstable. In what follows, exactly 
which type of modes (acoustic or collapse) do fall into the category defined by ()4.23j) will 
become apparent. 

From the results of static spacetime, we suspect that a Jeans instability must exist, 
and in such a region, it should follow that one of the Aj be zero, rendering ()4.23|) false. 
Consequently, we need to find the critical points of the Aj, dependent on the wavenumber 
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k. A consideration of the equation \\(x) = gives a solution for a "critical time" x — Xc- 

It turns out that A3 however has no such time. 

Let us examine the behavior of Xi(x) around x ~ Xc more closely. We set x 2 — Xc(l + e ) 
for a small parameter e, and expand Ai in powers of e. It turns out that 

/ + ^/FT^g = -l2eKlKl K Jf B l K Jf D + 0(e 2 ), (4.25) 

D e B + A B e D 

so that Ai oc \f--t for x ~ Xc- This dependence of Ai on e gives a clear picture of how A x 
changes around the critical point. For x > Xc, e > and Ai is imaginary. This corresponds 
to acoustic oscillations, in the stable part of A;-space. For x < Xc, £ < and Ai is real. 
This corresponds to an unstable part of /c-space, so that Xc is an indication of the transition 
through the Jeans instability. The time parameters may be defined so that initial time 
corresponds to ao = 1. This is because the explicit magnitude of ao is not determined by 
cosmology. By the definition x — a^ 1 / 2 , it is clear that x begins at 1 and decreases with 
increasing time. This gives us two ways of looking at the Jeans instability. One way is to 
consider the instability at a particular instant in time. For a particular time x, a subset of 
modes will be unstable for values of k for which Xc{k) > x ( we stress that Xc is a function 
of k). We may then consider what occurs as these modes evolve through time from this 
particular instant. The critical time Xc is fixed for any one mode, so that the modes which 
were originally acoustic will become unstable as x ~^ xt ■ Consequently more and more 
modes pass through the instability as the Universe evolves. The physical wavenumber k is 
of course dependent on time, thus the dependence of the instability on a time Xc shows the 
inextricable link between the wavenumber and time. 

We wish to relate these concepts back to the result discussed for static spacetime, and so 
must ask how the critical time Xc is related to the critical wavenumber ku of the mixture 
of components. In static spacetime we defined 

as the dimensionless parameter, indicating the relation of a mode to the instability at Km = 
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1. To place this quantity in an expanding Universe context, the substitutions 

W B - % v\k 2 - 6K%, 
X 



W D - ^, v D k 2 - QKl 
X 2 



are required. This gives Km the following form: 



Kll ~ e B /Kl + e D /Kl-ll- (427) 



It is explicitly seen here that the scale of instability changes with time, as was explained 
above. The analogy with the one-component case discussed in may be made here, 
where solutions were found in terms of the one-component Jeans wavenumber Kja^ 1 ^ 2 . 
With x — a -1 / 2 , we see that the quantity x^ 1 i n the two-component case is the exact 
analogy of Kj for the one-component case. 

Now that we have determined that Ai approaches zero in a particular region, it becomes 
clear that the WKB approximation will not be valid in this region x ~ Xc, from the condition 
(14.23)1 . Let us examine A 1 /Af in more detail, to determine its behavior over the whole of 
fc-space. An explicit evaluation using the definitions of / and g from ()4.6|) and ()4.7|) gives 

a; 3 



A? (Aa) 3 

where 



+ (4.28) 



Cd/?13 

In the ensuing discussion we use the quantity k to describe the comoving wavenumber 
dependence of the quantities involved. If we denote the numerator of 1)4.29)1 by N, it is simple 
to see that iV 2 < cd/3 2 3 for all k [see ()4.14j) ]. so that \F\ < 1. We also consider the limits of 
F for large and small k. For k — > 0, F — > 1 and for k — > oo, F — > — e B if K B > Kp., and 
F — > e B — to if Kb > K B . This means that asymptotically F is independent of k, and since 
\F\ < 1 for all k, A^/A 2 does not change sign. It now becomes clear that the magnitude of 
X'i/Xi is mainly dependent on the factor (Aix) -3 - We already know that as k — > ku-, Ai — >• 
so that in this region it has been confirmed that X'i/X 2 3> 1. Suitable regions where WKB 
might be valid must be sought far from the neighborhood k ~ ku- The dependence of 
on AiX leads us to suggest the WKB criterion in the amended form Aix 3> 1- This makes 
good physical sense, because if Ai is considered a frequency /inverse of a dynamical collapse 
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time, the WKB criteria requires that a large number of oscillations/significant change in Si 
occurs during an expansion time. 

A similar analysis needs to be performed for A 3 . We find 

|-(ASF (1 - n (4 - 30) 

This too demands the criterion A 3 x 3> 1. To find how A 3 behaves, consider firstly the k — > 
limit. This limit gives A 3 — > 0, so this region of /c-space is clearly inappropriate for WKB 
analysis. We note, however, that since A 3 ^ at k = k~M, this region should be checked 
more closely. After a careful examination of A 3 : 

A 3 X = V3^1-(K 2 B + K 2 D )x 2 -e D p 13 , (4.31) 

it becomes apparent that A 3 % 3> 1 only when both K B 3> 1 and K D 3> 1. This corresponds 
to the region k ^> k M , and so the region k ~ k M must be excluded from consideration as 
well. For A 3 , we are only left with the region k ^> ku as fulfilling the WKB criterion. For 
completeness, the same reasoning should also be applied to Ai. The k — > limit applied to 
Ai gives 

3 ' (4.32) 



(Aix) 3 2V6' 

which is not much less than 1, as is required to define it as a region amenable to WKB 
analysis. Thus the k ku region is inappropriate for Ai as well. 

In conclusion, k ^> kM is the only region for which the WKB approximation holds. We 
may summarize the methods available to analyze the two-component problem in various 
regions of A;-space by the following classification: 

k ^> ku The WKB method will give acoustic oscillations for all modes, with the rapidly 
varying part of the solution taking the form exp(±i|Aj|x). 

k k-M A Frobenius expansion (small parameter expansion) of the solutions needs to be 
developed. Some growing and decaying modes following a power law behavior may 
be expected for solutions corresponding to Ai, while some low frequency acoustic 
oscillations may be expected for solutions corresponding to A 3 . 

k ~ ku This region of parameter space is not accessible to analytic solution. Some numerics 
will be required to investigate this interesting region. 
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We will continue with the WKB analysis in this paper, and show how the equations ()4.13|) 
decouple. An investigation of the other regions of fc-space are taken up in 

We still need to check how f3\ and (3% vary, to ensure that all the are small. In particular, 
we need to consider whether relations of the form 

^«A„ i, j = 1,3 (4.33) 

Pl3 

are true. The analysis proceeds very similarly to that described above for the derivatives of 
Ai and A3, and a full description will not be given here. As a brief example, it can be shown 
that 

P>1 1 (4.34) 



A3A1 XiX 

for k ^> once again fulfilling the criterion Ai% ^> 1. Other cases follow similarly. 
Given that the relations ()4.33j) do hold, we can finally make the important statement that 
6i < Xj, Vz, j if k > k M - 

C. The Solutions 

Now that we have worked out the region of fc-space in which the WKB method produces 
valid leading order solutions, we proceed to derive these solutions by decoupling the equa- 
tions. To illustrate how the equations ()4.13|) decouple, we begin with an example. Taking 
the top row of the matrix equation, the following suggestive equation for fi may be written: 

/{ - (Ai - e x )f x = -£2/2 - e 3 / 3 - e 4 / 4 - (4.35) 

This may be treated as a first order inhomogeneous ordinary differential equation (ODE) 
for fi. The homogeneous portion has a simple solution 

"X 



fi ~ ci exp 



(Ai - e x )dx 

XO 



(4.36) 



with a constant of integration c\. This may be considered to be a first approximation to the 
solution, though it remains to be shown that it is the full leading order result. To evaluate 
the inhomogeneous portion of the solution of ()4.35|) . we require the first approximations for 
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the other fa as well. A corresponding analysis to that illustrated for fi yields 



f 2 ~ c 2 exp 
f 3 ~ c 3 exp 
fi ~ c 4 exp 



- / (Ai + ei)dx 



(A3 - e 7 )dx 



VII 



(A 3 + e 7 )d\ 



AO 



(4.37) 
(4.38) 
(4.39) 



When these are substituted into ()4.35|) . we achieve the rather complicated result 



h 



exp 



(Ai - ei)dx 



AO 



cn + C12 / cfxe 2 exp 



AO 



AO 



+ ci 3 / <ixe 3 exp 



AO 
X 



-Ai + A 3 + ei - e 7 )dx 



L"' XO 



1 / <ixe 4 exp / (-A a - A3 + e x - e 7 )dx )• 



(4.40) 



There are a number of integrals present here which need to be estimated to determine how 
the approximation is to proceed. 

As an example of a generic type of integral to evaluate, consider 



dx £2 exp 



-2 



\idx 



AO 



(4.41) 



Since the region of interest is k ^> fcjy, Ai is imaginary, and consequently the definition of 
the real function if)(x) = — *A(x) is useful. We integrate / by parts to obtain 



a; 



2 exp 



-2i / ^{x)dx 



AO 



AO 



AO 



d Aj 
Af 



exp 



-2i 



AO 



(4.42) 



It is our aim to show that all corrections to the leading order results (|4.36|) - (|4.39|) are of 
0(e), where e generically denotes any of the small quantities 6j, % — 1, 2, . . . 8. It is already 
obvious that the first term in (|4.42|) is of 0(e), because it consists of a \ x j term multiplied 
by a phase factor. We may constrain the second term by 

x 



AO 



d X [ ~± ) exp 



-2i 



^(x)dx 



AO 



< 



dx 



AO 



dx\\ 



Ai 
A? 



AO 



(4.43) 



This shows that indeed I ~ 0(e). 

The same type of analysis may be performed for the other more complicated integrals. 
In general, integration by parts is involved, and the terms may be shown to be bounded by 
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some expressions of the form X\/X 2 ^ or fi'JiP^Xj). It is found that in the WKB approxi- 
mation, ()4.36|) is the correct leading order term, and the integrals arising from considering 
the inhomogeneous portion of the ODE (|4.35jl are all of 0(e) below this leading order term. 
Similar correction integrals to those found in f)4.4()j) may be written down for f'2, fz and f^, 
all of which take the same generic form as those evaluated in the fi case. A lengthy analysis 
will show that all the corrections are of 0(e), so we have shown that to leading order we may 
effectively neglect the off-diagonal e« in (j4.13|) . In conclusion, the full leading order WKB 
solution to f is given by (|4.36)) - (|4.39|) . 

These solutions contain two important terms in the integrals. The Ai and A3 terms repre- 
sent the rapidly varying oscillatory portion of the solutions, i.e the control factor mentioned 
earlier. To obtain the true frequencies, these integrals need to be converted to integrals over 
t. We expect the e\ and €j terms to represent some slowly varying time dependent ampli- 
tude. To reveal the time dependent structure of the solutions more explicitly, the integrals 
of e\ and e-j need to be evaluated. This is facilitated greatly by the relation 



m = --■ 



For Oi consider the integral 
Using (|4.44jl once more we find 



01 



-dx = „ log 



xo Pi + cb/cd 
Pi(x) 



dx 



ilog(/3? + ^ 



Xo 



For notational expedience, we define the tilde quantities 



Pi(x)-PM 
? t 

Pi(x) 



"2 log 



Pi(Xc 



Pi(xo) -Ps(xo) 



1,3, 



PM) 
Hx) 

with the obvious property that f(xo) 
the integrals as follows: 



Pi(Xo) 
Pi(Xo) - As(Xo)' 

Hx) 



Xi(xo) 



1, 3, 



(4.44) 
(4.45) 

(4.46) 

(4.47) 
(4.48) 
(4.49) 



Xo 

x 

XO 



t\dx 
e 7 dx 



1, for any quantity fix)- This allows us to write 

~ Pl3 1 (4.50) 

(4.51) 
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Here the integral involving t-j was calculated using identical techniques as just illustrated 
for €\. 

At last the final form of leading order solution for f may be written down, with the 
explicit time dependent amplitude and rapidly varying oscillatory part: 



f\,2 ~ Ci >2 



f'3A ~ C34 



1/2 



13 



exp I ± / Aidx 

A3 



x 

exp ( ± / X 3 dx 
'xo 



(4.52) 



(4.53) 



We are now in a position to recover the original physical state vector x by multiplying these 
auxiliary functions by the matrix A, as given in the original definition ()4.4|) . The matrix A 
contains A's and /3's, which are not tilde quantities. The full general solution to ()2.6|) and 
fl2.7j) is finally given by 



Cl 



1/2 



1.3 



+ c 2 



+ c 3 



+ c 4 




Ciexp ( i / \\i\dx 
'xo 



| 2 exp ( -i I \Xi\dx 

xo 



| 3 exp ( i I \\ 3 \dx 

xo 



£ 4 exp 



Xo 



(4.54) 



Here the ^ are the eigenvectors defined by In summary, this leading order solution 

represents acoustic oscillations in the short wavelength limit, defined as 



k 2 > k 



4nGpB iiiGpD 



+ 



(4.55) 



B 



D 



which is a time dependent quantity. Equivalently, we may view the limit as given by 

% + % » (4.56) 

D 

A Jeans instability will not be evident for the solutions in this region of fc-space, but the time 
dependence will mean that the period of oscillation becomes longer until a point is reached 
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at which the WKB approximation is no longer accurate, and the solutions as displayed are 
not realistic representations of the underlying physics. Then different approximations need 
to be considered. The methods involved are discussed in detail in [11]. 



D. Relative Amplitudes of the Solutions 



The slowly varying time dependent amplitudes of the solutions 1)4.54)1 show how either 
dark matter or baryons dominate various modes, depending on whether HDM or CDM is 
being considered. This feature was noticed in the static case, and we now demonstrate it 
more fully through asymptotic analysis in the expanding universe scenario. The information 
is contained in the eigenvectors £ i( which give the relative amplitudes. It can be seen directly 
from (HOI that 

S f oc 1 = -f£ft (4.57) 
Od Pi 

for the Ai modes, and 

(4.58) 



5b 1 e D 

JT (X W = & 

Od P3 £b 



for the A3 modes. At first glance this may appear a little surprising, as there seems to be 
an asymmetry in the solutions. If the indices are interchanged D <-> B, the amplitudes do 
not appear to be the same, yet all such an interchange is doing is swapping the order the 
equations are written down in. 

As an illustration for the Ai modes 
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pi 




oc 




5b 




1 



with (3% 
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12e D 



(h + S), 



(4.59) 



whereas after the interchange D «-> B 



5b 




pi 




oc 




5d 




1 



with /3* 



X 



12e 



B 



-h + S). 



(4.60) 



The new quantity /3* is defined as being the new form of Pi after the interchange has been 
made. The square root term S, defined in (14.11)) is invariant under the interchange, whereas 
an examination of (|4.10[) shows that under the interchange h — > —h. It is however simple 
to show that (3i/3* = 1, which is to be expected from the symmetry of the differential 
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equations. The same result holds for the A 3 mode, where it can be shown that (3 3 (3^ = 1 for 
the quantities 



fo = £-(h-s), pi 



12e» 



-h-S). 



(4.61) 



We now examine the behavior of the amplitudes more carefully, to find the dominant 
components of matter. A useful large expansion parameter in the analysis will be the 
quantity y = \K B — K 2 D \x 2 - Let us also introduce the notation a = sign(K B — K 2 D ) = ±1, 
and expand the eigenvalues in y. For the Ai modes we may write 



Pi 
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2e D 



o + 



e D — e B 

y 



+ ^l + 2a(e D -e B )y- 1 + y 2 



(4.62) 



from which the following results may be deduced: 
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eD ~ e B 



+ o(r 2 ) 



, K B >K D {a=l) 



1 - (eg - eg) 

±e D y 



1 + 



eD — e B 

y 



+ 0(y 



-2\ 



K D >K B (a = -1) 



(4.63) 



This indicates that for CDM [K B > K D ) 



Pi oc (K B - K 2 D )x\ 

and the dark matter oscillations are dominant in the Ai modes, although as time increases 
they become less so. On the other hand for HDM (K D > K B ) 

and the baryon oscillations dominate, and their dominance increases with time. The A3 
modes display complementary behavior: 



•1 + (eD - e B ) 



03 



y_ 

eD 



^e D y 
1 - 



i + ^^ + oor 2 ) 



e B — en 



+ 0(y~ 2 ) 



, K B > K D (a = 1) 



, K D > K B (a = -l) 



(4.64) 



In this case for CDM 



P3 OC 



1 



{Kl-Kl) X ^ 

so that the baryon oscillations dominate, and the dominance increases with time, while for 
HDM 

P 3 oc -{Kl - K 2 B )x\ 
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which shows that dark matter oscillations dominate, but the dominance decreases with time. 
In summary, in a CDM scenario (Kp > Kb) baryons dominate the A3 mode and dark matter 
dominates the Ai mode, whereas the situation is opposite for HDM. 

The apparent asymmetry of these results with respect to interchange of subscripts B <-> D 
can be easily explained by taking into account the relations (|4.59j) - (j4.61|) . It was shown that 
after an interchange B <-> D, the new /3's are given by 

PI = -fo and f3l = -Pl (4.65) 

This is the correct way to view the asymptotic forms for the /3's after an interchange, rather 
than by a direct swapping of subscripts. By deriving the asymptotic results, information 
has been lost and direct swapping is no longer valid. As an example, for Kb > K^ 

fa ~ [K% - K 2 D ) X 2 - & = -ft ~ {K2 2 , (4.66) 

and the expected result P\fi\ = 1 holds. 

Now that we have gained some insight into the behavior of the modes in a generic sense, 
let us be a little more specific and investigate the full leading order solutions under some 
tighter physical constraints to see some more physical effects emerge. 



E. Initial Conditions and Resonances 



We attempt to reproduce the resonances described in the static spacetime case by im- 
posing some initial conditions on our solutions, and eliminating the arbitrary constants of 
integration q, % = 1, . . .4, using a similar procedure to that employed previously. In the 
static spacetime scenario we just stated the results, but here we explicitly go through the 
derivation. The general solution may be expediently written in the form 

x i(x) = ^2 c J v U)i(x) exp ( I Xjdxj • (4.67) 

j = l \Jxo / 

The amplitudes are given by the vectors 

/ ~ \ 1/2 / ~ \ 1/2 

v,1W2, = (|tJ ei * 2 ' V(3W4,= (llj 4l ' (468) 

At x — Xo, v (j)(xo) = £j(xo)> an d an equation for the initial conditions Xio is obtained: 

4 

x io = ^2 c £(j)i(Xo)- (4.69) 
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Note that henceforth any variable subscripted with a (possibly together with other sub- 
scripts) denotes that quantity evaluated at % = xo- As a reasonable simplifying assumption, 
we once again take x w = x 30 = 0, i.e. the perturbations start from rest. This gives a simple 
algebraic system for the q, with the solution 



(4.70) 



where the ratio of initial conditions Q = x i0 /x2o is used once more. 

When the above expressions for the Cj are substituted back into the general solution 
(|4.54|) . a simplified result follows: 
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|A 3 |dx 



(4.71) 



Let us concentrate in particular on the X2 and £4 components, which describe the actual 
matter content of the Universe, and may have some interesting implications for structure 
formation results. We may write the solutions in a form analogous to the static spacetime 
results: 



Sd(x) 
5b{x) 



X20 



£40 



Ci cos 
C 3 cos 



X 



AO 



\i\dx I + C2 cos 



\h\dx 



\i\dx I + C4COS 



AO 
X 



|A 3 |c?x 



Ad 



(4.72) 
(4.73) 
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The amplitudes are given by the expressions 



d = A^^(#V 3 , (4-74) 




P130 \ P3A3 

Cs = ^^f#V /2 , (4-76) 




C 4 = A °~ g ° 1 [^) . (4.77) 



P130 \ /3 3 A 
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The immediate obvious differences for these amplitudes with the static spacetime results 
are: 

1. Missing factors of ~ for d and (3, because the solutions are now all of a cosine form, 
rather than real exponentials (which is due to the fact that we are considering the 
k^> k M region). 

2. The amplitudes are all time varying. 

3. The amplitudes contain extra tilde factors which equal one at the initial time, but in 
general contain other time dependent terms not present even as constant factors in 
the static spacetime amplitudes. 

For x = Xo the amplitudes correspond exactly to those of the static spacetime amplitudes, 
and Figsl5HT2l are accurate representations of the amplitudes over a wide range of k. As 
time increases, the amplitudes tend to grow, but retain the same qualitative shape with the 
most marked features still occurring around the scale Km = Kmc- Since the WKB solutions 
are only valid for k ^> kM, this interesting region of fc-space does not apply, and the Q as 
presently defined should not be extrapolated to have any meaning around Km = Kmc- in 
the Km 3> 1 region the Q show no significant behavior, just tending to constant values close 
to zero or one. In particular no resonances are apparent. Thus the potentially interesting 

Q 

resonant features discovered by de Carvalho and Macedo |7J do not apply in the physically 



more realistic expanding universe scenario. 



from the small k solutions presented in 



11] 



Any significant effects would have to be sought 
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In addition to these statements, it must be added that the Aj eigenvalues characterizing 
the modes in the large k region do not have any physical significance in the small k region 
either. This was already apparent in the small k expansions for the one-component solutions 
discussed in The characteristic Jeans dispersion relation 
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V^P - AnGpo. (4.78) 



is not apparent in the small Kj expansions presented in [JJ|, and likewise, the A» found in 
the present paper are not apparent in the general solutions given by Eqs.(5.3)-(5.10) of jl]]. 
With these facts in mind, the discussion of the physical information we are able to extract 
from the WKB solutions at present is complete. 



V. CONCLUSIONS AND FURTHER WORK 



The structure and behavior of the eigenvalues and eigenvectors of two-component cos- 
mological density perturbations have been studied in great detail in this paper. We have 
reviewed the previous work done in a static spacetime background, and produced further 
results in this simple context. This has enabled the far more difficult expanding universe 
problem to be tackled. 

The WKB method employed has produced the full leading order behavior of all the modes 
in the Einstein-deSitter expanding Universe scenario. These solutions represent acoustic 
oscillations for wavelengths much smaller than the Jeans scale. The Jeans scale of the 
mixture has arisen in a natural way out of the analysis of the eigenvalues obtained through 
the WKB method, with some interesting interpretation. It is now a straightforward task 
to adapt the methods developed here to study a variety of further cosmological plasma 
modes. The ion-sound and two-component Langmuir oscillations would follow directly from 
the results presented here, and more complicated modes involving magnetic fields could also 
be obtained by similar procedures. 

We have also obtained the time- and fc-dependent amplitudes of the modes in a fairly 
general setting (the one restriction being initial perturbations beginning from rest). These 
results have shown that the amplitudes are very constant in the region of interest. The 
existence of resonances in the amplitudes found for static spacetime results do not apply 
here, as all resonances occurred for wavenumbers far smaller than ku- Thus a resonant 
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amplitude cannot be viewed as a mechanism for producing structures of a preferred scale in 
a two-component model. The eigenvalues derived in this paper do also not have any direct 
physical interpretation around the Jeans scale, or for small k expansions of the solutions 
of Eqs. (j2.6|) and (J2.7)) . Thus the results obtained in this paper must be considered to be 
restricted to the parameter regions considered here. 

It may be interesting to investigate different models such as a three-component HDM+- 
CDM+baryon fluid, or models involving a cosmological constant (especially given the weight 
of current observations Q|, [3], The analytics would become considerably more com- 

plicated, but some other interesting resonant scales may be found with a direct implication 
for structure formation. 
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FIG. 2: The eigenvalue Ai for a range of ed and cb, with V = 0.001 in the static universe scenario. 
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FIG. 4: The eigenvector amplitude functions (3\ and fa for HDM and CDM in the static universe 
scenario. 
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FIG. 5: The /c-dependent amplitudes of the Ai dark matter modes in a HDM universe for a range 
of initial conditions Qq = £40/2:20 • 
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FIG. 6: The /c-dependent amplitudes of the Ai dark matter modes in a CDM universe for a range 
of initial conditions Qo = £40/2:20 • 
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FIG. 7: The /c-dependent amplitudes of the A3 dark matter modes in a HDM universe for a range 
of initial conditions Qq = £40/2:20 • 
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FIG. 8: The /c-dependent amplitudes of the A3 dark matter modes in a CDM universe for a range 
of initial conditions Qo = £40/2:20 • 
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FIG. 9: The fc-dependent amplitudes of the Ai baryonic modes in a HDM universe for a range of 
initial conditions Qq = x<io/x2o- 
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FIG. 10: The /c-dependent amplitudes of the Ai baryonic modes in a CDM universe for a range of 
initial conditions Qq = X4q/x2o- 
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FIG. 11: The /c-dependent amplitudes of the A3 baryonic modes in a HDM universe for a range of 
initial conditions Qq = X4q/x2o- 
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FIG. 12: The /c-dependent amplitudes of the A3 baryonic modes in a CDM universe for a range of 
initial conditions Qq = X4q/x2o- 
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